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ABSTRACT 


A laboratory model of a single-axis vibratory angular rate servo has been 
developed, designed, and built. The vibrating element is a "free-free" metal 
beam mounted at its first mode nodal points and is vibrated by piezoelectric 
transducers mounted on the beam faces. Changes in the beam vibration pro- 
duced to Coriolis forces due to applied angular rates are detected by other 
piezoelectric transducers. Electronics are used to provide voltages to 
vibrate the beam and to amplify and demodulate the sensor voltage output. 

The major emphasis in the model development was placed on improving the 
threshold and null stability of the sensor. To achieve this, the new flexural 
pivot mounting method was developed which stabilized the drift characteristics 
of the sensor and improved the external vibration performance. Also, 
electronics were developed that provided a more accurate driving signal and 
a study was made of cross -couplings effects in the beam. The results of 
this work was at least at 5 to 1 improvement in sensor threshold and a cor- 
responding improvement in null stability. 
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FOREWORD 


This report was prepared by the General Electric Company, Avionic Controls 
Department, Binghamton, New York, on Contract Number NAS8-20618, 
"Microminiature Rate -of -Turn Sensor", administered by the National Aero- 
nautics and Space Administration, Marshall Space Flight Center. The con- 
tract was initiated under Control Number l-6-40-62363(lF) and S2(1F). 

This is the final report covering the sensor development work for the fifteen 
month period from June 1966 to September 1967. 

The contract work was a group effort, with the following principle contri- 
butors and their field of effort: S. B. Hamilton, Program Management and 
Analysis; J. N. Russell, Electronics Design; W. D. Gates and P. R. 
Danyluk, Mechanical Design; T. F. Conroy and G. Macko, Mechanical 
Details and Assembly. 

This report is published by the Avionic Controls Department under publica- 
tion number ACD 8569. 
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Section 1 


INTRODUCTION AND SUMMARY 


This report describes the development, design, fabrication, and testing of 
a laboratory model for a single-axis vibratory angular rate sensor. The 
vibrating element of the sensor is a "free-free" metal beam of approximately 
a square cross-section that is mounted at its first mode nodal points, and is 
vibrated by piezoelectric transducers. The transducers, bonded to each of 
the four faces of the beam, are also used to detect changes in the direction 
of beam vibration due to angular rates about the longitudinal axis of the beam 
which is the sensitive axis. Electronics, utilizing microelectronic opera- 
tional amplifiers, produce the voltages which excite the transducers, and 
also amplify and demodulate the sensor output signal. 

The impetus behind the development of vibratory angular rate sensors in 
general, and the "free-free" beam sensor in particular, is the potential 
long-life and excellent operating characteristics of this type of sensor when 
compared to conventional angular rate sensors with rotating elements. The 
absence of bearings and requirements for absolute dimensional stability, 
combine to make the expected life of a vibratory sensor much greater than 
its conventional counterpart. This, coupled with the accuracy, wide dynamic 
range, low power, dc operation, and wide bandwidth of the device, make the 
vibratory or "solid-state” sensor applicable for long duration space missions. 
The principle limitations of the "free-free" angular rate at the present time 
are its environmental performance, specifically temperature, and its 
"threshold" or ability to measure very slow rates -of-turn. The second of 
these was the subject of the work on this contract. 

A picture of the laboratory model is shown in Figure 1, and consists of the 
mechanical assembly of beam, mount, and transducers and the associated 
electronics. Table 1 is a summary of the model characteristics as deter- 
mined by test. The threshold/resolution of this particular unit is approxi- 
mately 0. 002 deg/sec and the null uncertainty is 0. 15 deg/sec. These 
parameters are an improvement by an order of magnitude as a result of the 
analysis and testing accomplished during this contract. The specific design 
changes to the "free-free" beam sensor which resulted in the improvement 
were the following: 

1. Development of the flexural pivot mount for the beam. This 

approach allows the interaction between the beam and mount to be 
minimized and held constant. 
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Figure 1. Solid-State Sensor Mechanical Assembly 
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Minimization of cross -coupling between input and output vibration 
to reduce sensor noise threshold. 


3. Development and design of improved drive electronics. The re- 
duction of harmonic content and improved regulation reduce sensor 
noise. 


Scale Factor: 

Full Scale Rate: 

Linearity: 

Frequency Response: 
Natural Frequency: 
Damping Ratio: 

Threshold and Resolution: 
Null Uncertainty: 

Output Voltage: 

Output Impedance: 

Power: 

Null Shift with Temp: 
(-10°C to +70°C) 

Null Shift with Vibration: 
lg (0 to 1500 Hz) 


0.395 volt/deg/sec 
±25 deg/sec 

< 1% (point-f or -point) 

1 db to 40 Hz 
63 Hz 

0.6 

~ 0. 002 deg/sec 
0. 15 deg/sec 
±10 volts dc 

< 500 ohms 

2 watts at ± 15 volts dc 

< ± 6 deg/sec 

< 0.25 deg/sec 


Laboratory Model Performance Parameters 
TABLE 1 
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Section 2 


THEORY OF OPERATION 

The operation of a vibratory angular rate sensor depends upon the motion of 
individual mass particles in a rotating frame of reference. The mass can 
be vibrating in any one of many modes (extensionally, bending, torsionally, 
etc.) at small amplitudes such that the particle motion is approximately 
linear and varying sinusoidally with time. Usually, the vibrational motion 
is constrained so that the mass particles can move only in one plane. The 
following paragraphs contain a description of particle motion and derivation 
of the vibratory sensor sensitivity. 

DERIVATION OF GENERAL EQUATIONS OF MOTION 

If the massjparticle position in a moving co-ordinate system with unit vec- 
tors, i, j, k is given by 


r = ix + j y + kz (1) 

arid if the motion of the i, j, k triad with respect to the inertial triad, I, J, 
K, is described by the angular velocity vector 

u> = ico + ico + k (2) 

x J y z ' 9 


then the velocity and acceleration of the mass particle are given by 


v 


dr 

dt 


3r 

dt 


+ co X r 


— dv 9 2 r n — „ 9 r Aw Xr + wXuXr 
a= dt = ^2 + 2c0X 9t + 9T 


(3) 

(4) 


3r 

where is the rate °f change of position with respect to the moving co- 
ordinate system. 


= ix(t) + jy(t) + kz(t) 


(5) 


4 


K the particle mass is m and the general force on the particle is F, then the 
equation of motion is 


1 — 9 r „ — 9 r 9 co — — — — 

i F =-2 t2 “ x »t t 5r X ""X" x r (6) 

3 1 


In Equation (6), the Coriolis acceleration is given by 

A C = 2 " x TT (?) 


Equation (6) can be expanded and re-written in matrix form as 



The specific case of Equation (8) that applies to the general vibratory angular 
rate sensor is the case when the particles are constrained to move in the 
y - z plane, and the angular rate u has only one component, w . Hence, 
the equations simplify to z 
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The generalized forces F and F consist primarily of internal damping and 

y ^ 

restraint forces of the vibrating material, and external forces produced by 
electrical, magnetic, or mechanical transducers. It is assumed here that 
the amplitudes of any movement are small enough so that the y and z damping 
and restraint forces are independent. Therefore, 


F 


f 


1 

* ^ 

u 

1 


K y 

y 

= 

y 


y 


y 

F 


f 


C z 


K z 

z 


z 


z 


z 


Transducer Damping Restraint 


Equations (9) and (10) combined yield the following particle equations of 
motion 



( 11 ) 


where 


C 

_y 

m 




K 

_JL = 
m 


it) 




( 12 ) 


The above equations are the basic operating relationships for a vibratory 
angular rate sensor. For convenience and purposes of definition, the *y’ 
axis is designated as the "vibrational" axis, the 'z' axis as the "output" axis, 
and the 'x' axis as the "input" axis. For example, if a free-free beam vi- 
brating in its fundamental bending mode is chosen as the vibrating mass, the 
sensor axis would be defined as shown in Figure 2. 
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Figure 2. Sensor Axes 
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The parameters for the vibratory sensor are chosen such that the centrifugal 
and Coriolis forces due to motion in the ’z* axis or output axis in the vibra- 
tional axis are negligible. 


Therefore, 


y(t) 

>> 

z(t) 

y(t) 

>> 

z(t) 

CM 

O 

3 

>> 

« 2 
X 

2 

OJ 

o 

» 

2 

00 
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The above conditions yield a lower limit for the possible mechanical resonant 
frequencies for a sensor with a given dynamic range. It will be shown in 
the following paragraphs that the sensor scale factor is proportional to 

GO 


— . If co is chosen too small, the output motion is significant for high 
c 


applied angular rates and the sensor output is non-linear due to the centri- 
fugal and Coriolis forces. For instance, mechanical resonant frequencies 
of at least 1000 Hz are needed for a sensor to be linear for angular rates of 
1000 deg/sec. 


SENSITIVITY TO A CONSTANT ANGULAR RATE 


If the moving co-ordinate frame is rotating at a constant angular rate co , 

z 


then the Coriolis acceleration of a particle is given by co y(t) where y(t) is 

z 


the amplitude of an externally induced vibration along the vibrational axis. 
The motion of the particle along the output axis is determined by this 
acceleration and the output axis dynamics of the particle. If the external 
force, f , is controlled such that 

y 


y(t) = y cos co t 
J ' J o c 


then z(t) is the solution to the equation 

. 2 

■ 2w z w c y o sin “c 1 = + Q~ + w o z ^ 
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The steady state solution to the preceding equation is 


z(t) = 


-2y co 

J 0 X 


CO 


CO 


2 . 2 


CO 


rrur sin (w c‘ - V 


1 - 


CO 


co Q 
c)V 


where 0 = arctan 
r c 


CO CO 
C 0 


~ , 2 2 , 
Q (co -co ) 
o v o c 7 


( 15 ) 


Equation (15) is the basic sensitivity expression for a vibratory angular rate 

sensor. The output motion is proportional to both the input or induced 

motion and the angular rate co , and is also dependent upon the relative 

z 

location of the two resonant frequencies, u> c and co q . 

BANDWIDTH - SENSOR RATE FREQUENCY RESPONSE 


The rate frequency response of a vibratory sensor with separated frequencies 

can be shown qualitatively and quantitatively to be approximately that of a 

lightly -damped second-order system with the natural frequency being the 

separation frequency, co . If the sensor is subjected to an oscillating input 

s 

rate of frequency, co^, about its sensitive axis, then the Coriolis force acting 

along the output axis is amplitude -modulated with a carrier frequency, co . 

c 

The force then has two sidebands, co + co. and co - co. which are amplified 

and phase-shifted by the dynamics of the structure in the output direction. 

If the sensor is designed for a frequency separation such that there is a 
minimal coupling (0 « 0), then the rate frequency response depends almost 

entirely on the upper sideband (co + co.) response. As the input frequency, 

co^ increases, the lower sideband is attenuated and remains at approximately 

zero degrees phase shift. As co c + co^ approaches co q , the upper sideband is 

amplified by the high Q of the output axis. The rate frequency response of 
the vibratory sensor therefore exhibits a resonance when co. = c o g due to this 

upper sideband amplification. Accordingly, the transfer function of the 
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where Q is related to Q . 
s o 

The equation of motion in the *z’ direction is given by 

oj „ f (t) 

z (t) + z(t) + 0 ) ^ z(t) = — — 
w Q v ' o w m 
o 

where f (t) is the Coriolis force, 
c 


In Laplace transform notation, Equation (17) becomes 


F (s) 

z(s) = G(s) - C » 
maj 


where G(s) = 



(17) 


If Q is large and the operating frequencies are near the resonant frequency, 





If the input rate and induced motion are 


CO (t) = CO cos co.t 
X XO 1 

(21) 

y(t) = y 0 Sin co c t, 

(22) 


then the Coriolis force f (t) is given by 

c 


f (t) = -2mx co co cos co.t cos co t 
c o xo c 1 c 


f (t) = -nw co cos (co + co.) t + cos (co - co.) t 

C J 0 XO C <- c 1 C 1 


For this force, the steady-state solution of Equation (17) is found to be 


(23) 

(24) 


z(t) = B u COS [( w c + CO.) t + 0 U J + cos [ (<c c - CO.) t + ^1 

where 

B u = • y ° 2 C B (“c + " t )> K = *<«> c + “j) (25) 

CO 

o 

B i = - y ° -- f C B (« c - COj), = 0 (o. c - co.) 


In order to provide sensor null stability as previously described, the reso- 
nant frequencies are separated such that co q > o> c and 0 (co c ) = 0. There- 
fore, the lower sideband phase-shift is always approximately zero. 

0. = 0 for all co. 
r £ i 
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Using this fact and expanding Equation (25), it becomes 


z(t) = B cos 0 cos co.t cos co t - B„ cos 0 sin co.t sin co t 
w u r u 1 c u r u 1 c 


-B sin 0 sin co.t sin co t - B sin 0 sin co.t cos co t 
u r u l c u r u i c 


+B. cos co.t cos co t + B„ sin co.t sin co t (26) 

£ l c S. l c ' 

A typical vibratory sensor uses a phase -sensitive demodulator to reduce 
harmonic noise and null shifts. The only components of z(t) that produce 
a rate output are those which are in quadrature with the induced motion y(t) 
(terms containing cos co t). If z (t) is the demodulated output of the sensor, 
then c 0 


s Q (t) = A(co.) cos [co.t - 6 (co.)] 


where 


A(co.) = (B cos 0 + B„ ) + (B sin 0 

v r v u r u £ ' v u r 


B sin 0 

<9(co.) = arc tan - = ^ — r — -^5- 

v 1 B cos 0 + B. 

u r u S. 


/] 


If the expressions for B u , B^ , etc. are substituted into Equation (27), then 
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1 


(31) 


6 (co^ = arctan 


2Q i ( ! 



where 



Q co 
o s 

CO 


(32) 


The gain and phase of the output given above are approximately those of a 
second order system with a resonant frequency of co and a quality factor Q. . 

S 1 

For example, if Q q = 750, co q = 1540 Hz, u> c = 1500 Hz, then co g = 40 Hz and 
Q. = 19.5. 

l 

DESCRIPTION OF SINGLE -AXIS FREE -FREE BEAM VIBRATORY 
ANGULAR RATE SENSOR 

GENERAL 


A single-axis sensor, consisting of a free-free metal beam mounted at its 
nodal points and vibrated at its fundamental bending frequency, is considered 
superior to other potential sensor configurations because of its basic simpli- 
city and its theoretical insensitivity to external vibrations or accelerations. 
Piezoelectric transducers, mounted on the four beam faces, are used to 
produce the motion along the vibrational axis and to detect the motion along 
the output axis. This transducer configuration minimizes the sensitivity of 
the sensor to dimensional stabilities between the beam and the external case. 
The following paragraphs contain an analytical description of the force 
transmission from the transducers to the beam, and the subsequent motions 
and voltages. 

The free -free beam and its attached transducers are depicted in Figure 3 
with the associated nomenclature. The driven motion along the vibrational 
axis is produced by applying a voltage to the drive transducer which exerts 
a force on the outside of the beam, which tends to contract or expand the 
skin of the beam between the nodal or mounting points. This force results 
in an effective torque on the beam about the nodes, and tends to bend the 
beam in its characteristic first mode shape as shown in Figure 4. If the 
frequency of the applied voltage is chosen to be the first natural frequency 
of the beam, then vibrational resonance occurs and the beam has the motion 
envelope as shown in Figure 4. If an angular rate is applied thru the mount- 
ing points to the longitudinal axis of the beam, the Coriolis forces produced 
by the velocity of the beam in the vibrational axis, cause the beam to vibrate 
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Figure 3. Free-free Beam & Transducer 





Figure 4. Beam Vibration Envelope 
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along the output axis. This vibration bends the output piezoelectric trans- 
cuder which produces an ac output voltage whose amplitude is proportional 
to the applied angular rate. Typical vibration envelopes for the free-free 
beam are shown in Figure 4 for both the vibrational axis and the output axis. 

INDUCED BEAM VIBRATION ALONG VIBRATIONAL AXIS - 
TRANSDUCER FORCE -VOLTAGE RELATIONSHIPS 

The induced beam vibration as referred to in the previous sections is pro- 
duced by applying an ac voltage to the "drive” piezoelectric transducer. The 
effective force produced between the ends of the transducer is given by 


F 


DR 


W 


K e g 31 


DR 


(33) 


where W is the width of the transducer, 


V 


DR is the applied voltage, is 


the transducer 'g' constant, and K is the transmission factor which takes 

6 

into account the mounting method for the transducer. This force produces a 
torque on the beam about the mounting points because of the beam thickness. 
Referring to the nomenclature of Figure 3, the torque is given by 


, _ 2 4 Fh 

DR 3 L 


(34) 


The angular deflection of the beam about its mounting points away from its 
equilibrium position can be approximately described by the differential 
equation 


— 1^-. 6 + C 9 + Kt-, 0 
Y B y 


DR B y 


(35) 


where I R is the inertia of one -half the beam about the mounting point and 
Kg is the rotational spring constant of the beam. Therefore, 


e (s) 
y 


t dr^ 

k b 


i + 


co Q 

y y 



00 


y 


(36) 
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where 


co 


2 


y 


V 


co Q 

y y 



(37) 


The deflection angle 9 is shown in Figure 4. The quality factor is usually 
very high (~1000) and represents only the material damping effects. 

If the deflection angle, 0^(t) is known, then both the maximum deflection of 

the beam at the center and the extension of the "feedback transducer" can 
be found. These are both important parameters since the first determines 
the sensor scale factor and the second influences the drive circuit design. 

If the first mode deflection curve of the beam is assumed to be approximately 


y(x, t) = sin 



(38) 

(See Appendix A) 


then the mode deflection angle 6 


y 


is related to y Q by the equation 


o 


h*k e 

TT y 


and 


y o (s) = 


2.6 1 h w „ 

3 77 K fi g 3 i V DR 


, s s 
1 + — — + 


(s) 


co Q 

y y 


CO 


(39) 


(40) 


If the drive voltage is designed to be generated at the resonant frequency 
co , Q and 

y y 


V DR (t > " V DRO C0S V 


(41) 


then 


y 


o 


2.6 1 h w Q 

31 Vsi y 


(42) 
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Equation (42) gives the maximum beam drive deflection in terms of the drive 
voltage amplitude applied to the drive transducer. When the beam is driven 
along the vibrational axis using the drive transducer, then the "feedback" 
transducer (transducer on the opposite side of the beam from the drive 
transducer) is extended and contracted as the beam moves. This extension/ 
contraction produces a voltage on the feedback transducer that is proportional 
to the amplitude motion, and may be used electronically for amplitude and 
frequency control of the beam motion. The feedback voltage is related to 
angular deflection by 


V M (t) = 


th 


FB Ld 


31 


9 (t) (refer to Figure 3 for nomenclature) 


and therefore. 


2H th w 


3L K K d 31 g 31 

V FB^ - 2 V DR^ 


1 + 


+ 


u) Q 2 

y y » 


(43) 


(44) 


This is the steady-state transfer function on "gain" of the beam-transducers 
along the vibrational axis. It plays an important part in the evaluation of a 
given sensor and the design of the electronics that produce the drive voltage. 

OUTPUT AXIS MOTIONS DUE TO ANGULAR RATE 


With the above information, the beam deflection along the output axis for a 
given applied angular rate may be calculated as well as the corresponding 
voltage output of the readout transducer. For a constant input angular rate, 
£2 , about the x-axis, the Coriolis force acting on the beam per unit length 

along the beam is given by 


f cz (t) = 2|i« x y(t) 


(45) 


where 


y(t) = u) y y Q (sin ^ - |) cos u> y (t) (46) 

or 

f ct(t) = 2p£yo y y o (sin ^ - |) cos o> y t = A(t) (sin ^ - |) (47) 
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The corresponding torque in the output plane about the mounting points or 
nodes is found by integrating this unit Coriolis force. 

L/2 L/2 

T z =J x f cz (x, t) dx = A(t) /x(sin^-i) dx (48) 

0 0 


This results in 


T z (‘) 




ixfi co y 
x y J o 


cos co t 

y 


(49) 


The angular deflection response in the output axis is described by a differential 
equation similar to the one used for the vibrational axis. 


T (t) = I D d (t) + C6 (t) + K n 0 (t) 
z B z z B z 


or alternately 


«,w ■ 


K T - '■> 


- s s 
1 + — — + 


U) Q 2 

Z Z CO 

z 


(50) 


(51) 


where 


co 

z 


2 


K 


B 


B 


co Q 
z z 



If only the steady-state (or sinusoidal) operation is considered and the drive 
frequency co is significantly different from the output axis resonance co , 

y z 

then the peak angular deflection and the torque are given by the approximate 
relation. 



( 52 ) 
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The maximum beam deflection is related to angular nodal deflection by 


z 

O 77 


9 


if the first mode beam deflection is given by 

sin 7 tx 2 


Z ( X , t) = 

Therefore, 


J o L 


77 


COS (JI) t 

y 


z 

o“ 3 

77 


3 r 


-i] 




K 


B 


/ CO \ 3 


where y Q is related to the drive voltage 


( 53 ) 


(54) 


(55) 


The bending of the beam along the output axis also bends the "readout trans- 
ducer" and produces a voltage proportional to the amount of bending, and 
thus, to applied angular rate. The readout voltage V^q is related to the 

angular deflection in the readout plane by 


th 9 

V RO (t) = L K d„ 1 
e 31 


(56) 


from which the sensor scale factor (volts/deg/sec) may be found. 

The relationships established in the previous paragraphs serve the purpose 
of a foundation from which to study sensor operation from the point of view 
of developing design criteria and analyzing possible sources of sensor null 
shift, noise, and drift. 

The important facts that should be remembered from the preceding text are: 

1. The main component of sensor vibratory motion is at mechanical 
resonance. 

2. The output or angular rate sensitive motion is of resonance and 
approximately in quadrature time -phase with the main component. 

The following paragraphs examine the causes of sensor null shifts and noise, 
in general, and then examine the electrical and mechanical design of the sen- 
sor in detail. 
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SOURCES OF SENSOR NULL SHIFTS AND NOISE 


The primary sources of vibratory sensor null shifts and noise are directly 
related to or can be represented by changes or anomalies in the piezoelectric 
transducers or external distrubances such as mount vibration. For instance, 
if the force or torque produced by the drive transducer is not exactly aligned 
with the vibrational axis, then beam vibration will occur along the output 
axis that cannot be distinguished from the vibration produced by an applied 
angular rate. Factors such as uneven bonding between transducer and beam 
or a nonhomogeneous transducer can cause the above. The major sources 
of sensor drift and noise are listed below. 

1. Variation in the direction of drive transducer force with time or 
environment. 

2. Quadrature coupling between the vibrational and output axes. 

3. Beam vibrations induced by external accelerations. 


The sensitivity of the sensor voltage output to changes in the direction of 
drive torque or force may be represented by differences in applied voltage 
for parts of the drive transducer or differences in coupling factors. The 
force -voltage relationship 


F 


DR 


W V 


DR 


Kg 


e & 31 


(57) 


was used in the previous work in this section as the basis for sensor 
operation. If the coupling factor K e is not constant across the width of the 

transducer, then there is torque produced in the z (output) axis as well as 
the y (vibrational) axis. For simplicity, this situation may be represented 
by a split transducer with the two halves having different applied voltages 
as depicted in Figure 5. 

The total torque acting on the beam in the output axis is 

T z= T 1- T 2=36^[ V 1- V 2] < 5 ' 


and the resulting output voltage from the readout transducer due to beam 
bending is 


th T 


RO 


2L Wsi 



(59) 
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( 60 ) 


th W 


L V 1 - V 2 ] 


RO 


72 L K e K B g 31 d 31 




This output voltage is either in phase or out-of -phase with the drive voltage 
and therefore, cannot be distinguished from an output voltage due to an input 
angular rate. Therein lies the sensor null shift problem. Changes in trans- 
ducer parameters such as coupling factors, d and g constants with environ- 
mental and time, may be represented by changes in the effective applied 
voltage. The actual magnitude of this sensitivity, which will be calculated 
in the Sensor Design, Section 3, in connection with a self-test method, is 
large enough to pinpoint this sensitivity as the major limitation of sensor 
performance. Sensor parameters such as threshold and null stability are 
affected by this drive transducer stability problem. 


The quadrature coupling, between the vibrational and output axes, makes the 
sensor null voltage output sensitive to changes in the amplitude of the driven 
vibration. As was noted in a previous discussion, the output motion is ap- 
proximately but not exactly in quadrature with the input motion. The degree 
of exactness of this phase relationship depends upon the difference of the two 
natural frequencies and the associated quality factors. From the sensor gain 
equation given in this Section, the output motion of the beam drive to an angu- 
lar rate can be approximated by 


z(x, t) 


y ^ 
J o x 

CO 



2 

7 T 



( 61 ) 


where y Q is the amplitude of the beam motion in the vibrational axis, f2 x , 

is the input angular rate, and u> if the beam separation frequency. 

s 


The principal causes of crosscoupling between input and output are the small 
phase shift indicated in the above expression and the actual response of the 
output transducer. The transducer that senses the output motion cannot be 
made completely insensitive to the motion in the vibrational axis. Because 
of the sizeable difference between the amplitudes of input and output motion 

4 

(10 ) this sensitivity can be extremely important when considering sensor 
performance at null or very low rates. The output of the readout transducer 
may be represented as 


V c (t) = K z z(t) + K y y(t) 


(62) 
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where K z is the transducer scale factor discussed previously, K is the 

quadrature coupling coefficient, and y(t) is the beam motion along the vibra- 
tional axis. Therefore, if 

y(t) = y Q cos co t, (63) 


then 
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If the phase shift represented by the term „ ^ is small, so that sin _ ^ 

y co y co 
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= ^ y ■ , then referencing to the phase of the output signal, 
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where 


t, = t - 
1 Q co 


z s 


( 66 ) 


The resonant mode coupling and the potential cause of null shift or instability 
is the dependence of the amplitude of the sin co^t, term on y Q . If the null 

shift of the sensor due to this coupling is to be less than an equivalent angular 

rate, , then 
x 


y ^ 
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Null shift due to 

coupling (68) 
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Note that the preceding relationship is independent of the separation fre- 
quency co , as long as the approximations are satisfied. These approxima- 
s 

tions state that the phase shift of the beam motion in the output axis with 
respect to the drive voltage is small, typically less than 5 degrees. The 
magnitude of the coupling null shift can be minimized by either trimming the 
readout transducer or by splitting the readout transducer electrode and null- 
ing the voltage. However, even though trim is possible, the basic coupling 
still exists and changes in the drive motion amplitude are likely to be reflected 
in sensor null shifts. The same is true for drive motion amplitude noise 
which has bearing on the sensor threshold or S/N ratio if coupling is present. 
The optimum setting for the separation frequency is usually about 20 to 30 Hz 
for a beam natural frequency of 2000 Hz and a Q of 1000. 

The sensitivity of the beam motion to external accelerations introduced 
through the mounting structure. This sensitivity represents coupling that 
is a departure from true free-free beam vibration. One of the primary 
reasons that the use of a free-free beam as a vibratory angular rate sensor 
should realize better performance than previous vibratory sensor designs is 
the theoreticle insensitivity of the primary beam motion to accelerations 
introduced from the mount. As a result of the beam mode shape, as shown 
in Figure 6 , and the fact that the modes of the first mode are used as 
mounting points, the instantaneous linear momentum of the beam is zero at 
any time ( ends going down, middle up). If perfect nodal mounting is achieved, 
the stress on the mounts due to the transducer -driven beam motion is zero, 
and the interaction between the beam and the external environment is zero. 
The following analysis is intended to determine the mounting tolerances 
necessary to minimize vibration sensitivity to explore possible effects on 
sensor null stability. 

Referring to the nomenclature of Figure 6 , the neutral axis of the beam 
motion is represented by 

y(x, t) = £y o sin -a J cos co y t (69) 

where cu is the first free-free mode bending frequency. The nominal loca- 

y 2y o 

tion of the modes of vibration are at x = e =0. 220£ and a = a = — — . At 

O 0 7 T 

this point, the configuration is balanced and the net vertical momentum is 
zero. 

a 8 

M(t) = p. f qj- dx = 0 where p is the mass per unit (70) 

0 length of the beam. 
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If the beam is mounted at other points such as 


x = e = e + Ae 
o 


where a = a + Aa 
o 


( 71 ) 


then there is a resulting unbalance momentum. The modified motion and 
the resulting momentum are given below 


y(x, t) = y 
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M(t) = 77 py co cos — j?- Ae sin co t 

® y y 


(72) 


The total acceleration on the supports due to this momentum is 

“T = h y o cos ~r (t) cos “y l < 73 > 


For example, if co y = 2000 cps, y Q = 0.005”, £ = 4.0” 

and the mounting points are each 10 milli-inches (Ae) inside the nominal 
positions, then the total peaks acceleration on the mounts is approximately 
12.0 g's. Also, for the example the peak acceleration of the beam at its 
center is approximately 750.0 g’s. 

The representative numbers given above point up three potential problem 
areas with regard to the fabrication and assembly of the vibratory sensor. 
First, if perfect nodal mounting is not achieved, the sensor will be sensitive 
to external vibration at its resonance frequency (about 2. 0 deg/sec/g of 
2000 cps vibration in the readout plane for the above example) and, probably 
more important, stresses are set up in the mounting structure which may 
effectively shift the plane of driven vibration and result in significant null 
shifts of a random nature. Also, the above serves to accentuate the inade- 
quacy of a rubber mount when a high performance sensor is required. 
Second, at the beam accelerations that are necessary for operation, the 
stability of the transducer bend is in question. Any change in uniformity or 
strength will result in large permanent null shifts. Third, if the beam is 
not balanced, the neutral axis of vibration for the drive and readout planes 
will not coincide and the nodes for each plane will be in different locations. 
This fact would complicate the mounting problem. 
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COMPLETE MATHEMATICAL MODEL FOR THE FREE -FREE BEAM 
ANGULAR RATE SENSOR 

All the sensors parameters and considerations can be combined into a single 
mathematical model for the free -free beam vibratory angular rate sensor. 
The ultimate objective of the model is to describe the sensor accurately and 
in doing so predict the limitations on sensor performance. The results ob- 
tained from the hardware during tests will aid this prelimi- 
nary effort. The following discussion is intended to be qualitative and to 
define the various aspects of the problem. 

The significant segments of a mathematical model for the free -free beam 
sensor necessarily contain the following: 

1. Two dimensional beam dynamics described in fixed beam 
coordinates. 

2. Force-voltage-current relationships of piezoelectric transducers 
described in transducer coordinates . 

3. Electronic transfer functions and noise. 

4. Coordinate transformations from transducers to beam as a 
function of environment. 

5. Effects of mount on beam dynamics. 

6. Transmission of external accelerations through the mount to 
the beam. 

A block diagram of the mathematical model is shown in Figure 7. The "K" 
constants’ represent the transducer gains with associated geometry; the 
"C’s" are transducer coupling constants that are due to misalignments and 
anomalies; the "H’s" are the transfer functions of the electronic circuits. 

It is evident from Figure 7 that sensor performance is directly related to 
the stability and magnitude of the coupling C’s since they represent the only 
interactions between drive and readout beam motion besides input angular 
rate 0. The beam transfer functions X(s) and Y(s) are limited to first mode 
dynamics and it will be shown in later work that this representation is valid. 
External distrubances that act upon the sensor are represented by the noise 
inputs N and the force inputs F (external accelerations). 
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The various parameters of the model can be evaluated by one or more of 
the following tests: 

1. Open-loop tests of beam and transducers with an audio oscillator 
(X, Y, K). 

2. Calibration of electronic circuits (H, N). 

3. Static tests of sensor and electronics - transient responses or 
modulation (C). 

4. Performance tests of sensor in a controlled (K) environment. 

5. Environment evaluation of sensor F, • 

6. Analysis to extend experimental results. 

All of the transfer functions and constants have been developed in the previous 
paragraphs. 
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Section 3 

SENSOR DESIGN AND ANALYSIS 


The following paragraphs contain a description of the various design para- 
meters of the free- free beam vibratory angular rate sensor and the criteria 
by which they are selected. The discussion is broken into mechanical and 
electrical sections and is intended as the background for the specific descrip- 
tion of the laboratory model and its performance contained in Section 4 . 

MECHANICAL ANALYSIS AND DESIGN 

A sketch of the mechanical design for the laboratory model is shown in Figure 
8. The mechanical part of the sensor consists of the beam, the mounting 
structure for the beam, the four piezoelectric transducers, the bond between 
the transducers and the beam, and the electrical leads connected to the trans- 
ducers and brought out to the electronics. The design of the mechanical part 
of the sensor is a combination of material and parameters selection. 


One of the most basic selections for the sensor is that of the first free-free 
bending mode frequency of the beam and thus, the sensor operating frequency. 
It is shown in Appendix A that this resonant frequency is given by 


_ 22.4 / El \ 1/2 
- L 2 W ' 


(74) 


where El is the flexural rigidity of the beam, i u is the mass of the beam per 
unit length, and L is the length of the beam. The cross-sectional moment of 
inertia 0 , I, is given by 


I 



(75) 


if the beam cross section is square and of side h. Therefore, 


(76) 


and entailed in the frequency selection is the selection of the beam shape 
(length -to -width ratio). 
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Figure 8. Modified Gimballed Beam and Mount 
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The general considerations for this selection are the following: 

1) Sensor scale factor: Scale factor decreases approximately as 

2) Bandwidth of expected external disturbances: The resonant frequency 
should be above the expected vibration environment 

The selection of beam length- to- width ratio for a specific frequency is depen- 
dent on scale factor versus mounting considerations. If the cross section 
is increased, the sensor scale factor is increased because of the greater 
transducer bending; however, and increased cross section makes it more dif- 
ficult to obtain accurate nodal mounting because of the greater surface area. 


The separation frequency, or the difference between the vibrational and output 
resonant frequencies, is determined by making the beam cross section rec- 
tangular instead of square. This alters the cross-sectional moment of inertia 
(or the two axes). Consequently, 



Equations (77) and (78) may be combined to yield an expression for the ratio 
of the two resonant frequencies. 



(79) 


The tolerances on the beam cross-sectional dimensions are established by 
selecting a tolerance for the sensor separation frequency and then, deter- 
mining the dimensional tolerance needed to maintain this. 

The analysis presented in Section 2 indicated that the sensitivity of the sensor 
to external disturbances was dependent upon how close the beam mounting 
points were to the actual nodes of the first free-free beam mode of vibration. 

In fact, one of the reasons that the free-free beam is considered the best 
mechanization for a vibratory angular rate sensor is its theoretical intensitivity 
to external accelerations. When mounted exactly at the first mode nodes, 
the amplitude of beam vibration at resonance does not depend on mount accel- 
eration because the net beam momentum is always exactly zero. For this 
reason, second-order effects that result in nodal shifts are of prime interest 
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in sensor design. The mass of the piezoelectric transducers which are bonded 
to the beam are one such effect. The nodes of vibration are drawn in signifi- 
cantly by this added mass at beam center. Referring to Figure 9, which 
shows schematically the beam and transducer, the transducers are centered 
on the beam and are of length i . If the general first mode shape is given 
approximately by 


Y(x) = Y sin ~ - a (80) 

v ' o L 

Where "a" is the constant which defines the nodes of vibration, then the nodes 
for a beam without transducers are at 


x 


n 


and 


x 


n 


7 sl ” '*(#) s °- 22L 



0.78L 


(81) 

(82) 


If transducers are added, the beam is no longer uniform and the extra mass 
tends to pull the nodes toward the center. Since the nodes are defined as the 
mounting points for which the net momentum is zero, then Momentum of 
Beam + Momentum of Transducers = 0, which can be expressed approximately 
by the equation: 

L+4 

L 2 (83) 

M B f Q (y 0 sin -X " a ) dx + 4M T fh-L (y o Sin ¥ " a ) dx = 0 

~2 

Carrying out the integration and solving for the constant "a", 


M B + 4 M t 
M s + 4 (|;) M T _ 

where Mg is the mass of the beam 

M^ is the mass of one transducer 
l is the length of the transducer 
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(b) Gimbal With Dimensions 



Figure 9. (c) Beam With Transducers 
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The node location is found to be: 


x 


* 


n 


0.6366(M b + 4 M t ) 

M B + 4 (r) M T 


(85| 


As a general rule, the node shift due to transducer mass is approximately 
2 to 3% and is particularly important when the beam is hard-mounted such as 
the gimballed beam described in the Control Loop and Electrical Analysis 
and Design test that follows. 

CONTROL LOOP AND ELECTRICAL ANALYSIS AND DESIGN 


The electrical analysis and design for the free- free beam sensor consists of 
the definition of the control circuits needed to generate the drive voltage, to 
process the voltage on the readout transducer, and to provide an electrical 
damping signal for the beam. The gains, frequency responses, and non- 
linearities of the electronics are directly dependent on the mechanical beam- 
transducer parameters such as resonant frequencies, and are chosen to 
minimize sensor output noise and null shift, as well as to fix sensor perfor- 
mance parameters such as scale factor and damping ratio. The electronics 
are, of necessity, analog because of the amplitude- modulated beam output 
and are synthesized using an operational amplifier mechanization. 

DRIVE LOOP ANALYSIS AND DESIGN 

The drive electronics, or drive loop, which provides the drive voltage to 
vibrate the beam in the vibrational axis , can be generally designed by one of 
two methods. The first method, which consists of a fixed amplitude- variable 
frequency oscillator, uses the output of the feedback transducer to determine 
the frequency of oscillation and to adjust the oscillator frequency to the beam 
resonance. The second method consists of an analog closed loop oscillator 
for which the beam-transducer combination is the frequency determining 
element. As was noted previously, the so-called "gain" or transfer function 
of the beam-transducers with the drive transducer as input and the feedback 
transducer for output is very similar to a highly resonant series R-L-C circuit. 

For most vibratory sensor designs, the second of the above methods is the 
most desirable because of the tolerances and variations of the beam-transducer 
gains and resonant frequencies. Because of the self-tracking nature of the 
closed- loop oscillator, variations in sensor performance due to aging and 
environment are minimized by this method. If the resonant frequency of the 
beam-transducers should shift the oscillator frequency would track it because 
it is the frequency-determining element. Gain variations can also be com- 
pensated for, as will be seen in later paragraphs, by providing a loop gain 
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much greater than unity, and designing a saturating nonlinearity or regulator 
control to fix the amplitude of oscillation independent of the beam . The fol- 
lowing paragraphs contain a description of the drive loop dynamics as well as 
the effects of beam -transducer variations such as resonant frequency and 
gain, and the effects of drive loop noise. 

A block diagram of the basic oscillating drive loop is shown in Figure 10 and 
consists of the beam and transducer dynamics, and the drive electronics. 

The two functions provided by the drive electronics are filtering of the feed- 
back transducer voltage to assure the correct loop phase for oscillation, and 
an amplifier to assure a loop gain greater than unity under all conditions. 

The filter is chosen to be a double pole at the beam resonant frequency instead 
of an integrator to insure that the loop will not oscillate at any of the higher 
beam modes (Appendix A). The frequency stability of this loop in terms of 
its ability to track the beam resonant frequency is dependent on the difference 
between the beam resonant frequency and the frequency at which the filter 
provides -90° phase shift. The phase of the beam- transducer transfer function 
is 


0y(^>) = n - tan 


CO CO 

y 


[_Q (cd 2 - CO 2 ) 

L y y 


( 86 ) 


and the filter pahse shift is 

0 f (w) = - 2 tan -1 (87) 

where co^ is usually chosen to coincide with co^. The slopes of these phase 
functions at their respective characteristic frequencies are 
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Figure 10. Basic Oscillating Loop 
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If w, and are close together then both phase characteristics may be approx- 
imated by straight lines in this region 

*y <"> 3 % < W y) + Vdi) o> = Uy ( " ' “V 

*f<“> 3 + (S) (W '"f ) 

7T 7T 

Since 0 (co ) = ^ , and 0, (co^ ) = - -s- and the fact that the loop shown in 

J J 

Figure 8 oscillates at a frequency co where co is determined from the 

OS OS 

equation 

*y ("os* + *1 ( “os ) = 0 ( 9 °> 

provided that the loop gain is greater than unity. Substituting and solving for 

co , we find that 
os’ 


for high Q's (91) 

1 + y 

1 2Q co r 
T f 

If the quality factor is approximately 1000, then the frequency of oscillation 
is always within 0.004% of the natural frequency of the beam, if the difference 
between the natural frequency and the filter frequency is less than 5%. 

The loop in Figure 10 not only sets controls for the frequency of operation, but 
it also determines the start-up time for the sensor. A good rule of thumb is 
that for effective loops gains greater than 10, the start-up time is controlled 
by the beam Q (bandwidth). The amplitude of beam movement rises with an 
approximate time constant of t 
2Q 

r s y 

S O) 

y 
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The previous paragraph describing the oscillator operation did not consider 
the amplitude of vibration. For a loop such as that shown in Figure 10, the 
amplitude is determined by the saturation characteristics of the electronics. 
This saturation produces higher harmonics in the drive voltage that tend to 
distort the sensor angular rate output and raise its noise content. Also, if 
saturation characteristics are used for amplitude control, the sensor scale 
factor which is proportional to the amplitude of beam motion in the vibrational 
axis becomes very sensitive to changes in the drive transducer and bond. In 
order to eliminate these problem areas, the drive loop is designed with an 
automatic gain control which fixes the oscillation at a precise amplitude and 
maintains sinusoidal operation. A block diagram of this amplitude control is 
shown in Figure 11, and consists of the basic frequency-determining loop plus 
an automatic gain control for the drive amplifier . The amplifier gain is con- 
trolled with a signal proportional to the difference between the amplitude of 
the feedback voltage and a chosen reference voltate (e). The gain variation is 
specified as 

K(e) = (me + 1) (92) 


When m is the gain slope and Kg Q is the nominal value of beam gain. Using 
the above and the nomenclature of Figure 11, we find that 


V D (t> 


[l-“H<V R -Vo>] V> 


(93) 


where and are the amplitudes of the drive feedback voltages respectively, 
and V q is the reference voltage. 


Since it is much easier and more accurate to establish a dc reference instead 
of ac, the functional block H is of necessity a peak detector and low frequency 
filter, and may be represented by x 1 

H(s) = TTT? (94) 

Since the control loop is nonlinear, stability must be considered in the limits 
and particularly about the nominal operating point. If 
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Figure 11. Amplitude Control Block Diagram 
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then 


(97) 


where 
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Then, the total electronic transfer function for small changes in oscillation 
amplitude is 


av d , . JL 

- V R <S) ' K B 
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(100) 


The beam-transducer amplitude transfer function consists of a single log at 
the bandwidth of the resonant peak. Amplitude modulation produces sidebands 
around the resonant carrier which are attenutated the difference between the 
sideband and the carrier becomes greater. Therefore, the loop of Figure 11 
may be written as 


A(s) 


(m V Ro - 1) 


m Vro - 1 


U + T g ) (1 + T 0 S) 


( 101 ) 


A Bode diagram for a typical loop is shown in Figure 12. The RHP zero that 
is introduced by the filter and variable gain amplifier is always a higher fre- 
quency than the filter pole, and selecting parameters for adequate stability 
is not a problem . The regulating capability of this loop are dependent upon 
the dc loop gain (m - 1 ) . The actual regulation in percent in 

— — _ ; for a 30% change in beam gain and a loop gain of 10, the feedback 
mv Ro 

voltage amplitude change is approximately 3%. 

One of the primary design and analysis considerations for the drive oscillator 
control loop is its susceptibility to noise, either electrical or mechanical, 
that could be transferred output axis and affect sensor threshold/resolution. 
As was shown in Section 1 on Theory of Operation, the drive transducer inter- 
face is very sensitive to voltage changes that are reflected in sensor null 
shifts. Also, the possibility of some noise content at the sensor natural fre- 
quency puts a premimum on noise considerations. The noise response of the 
loop shown in Figure 11 with the noise simulated at the gain control input of 
the AGC amplifier is 
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Figure 12. Variation Amplitude Control Bode Diagram 
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If the time constant of the filter H is chosen to yield a critically-damped 
characteristic, then: 


av d (s) = 
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and 
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The loop gains and noise responses for various values of m E^. q are shown 

in Figures 13 and 14 for Tg =0.08 sec and Kg = 1 . These curves point up 

the fact that the drive loop is susceptible to broad band noise . Random mod- 
ulations of the drive voltage are especially significant if they excite the read- 
out plane resonance (noise at the separation frequency). 

The various typical curves given show that for high loop gains (good feedback 
voltage and scale factor regulation) that the noise response of the loop degrades 
especially at higher frequency. This indicates that there is a tradeoff and an 
optimum setting for regulation versus sensor noise and threshold. This fact 
will be borne out by the cross coupling data presented in the Testing Section 
when the sensor is subjected to controlled noise in the drive loop. 

READOUT AND DAMPING LOOP ANALYSIS AND DESIGN 


The readout or output functions of the sensor consist of amplification and de- 
modulation of the beam- transducer voltage output due to an applied angular 
rate . Since the beam-transducer system is an ac system with the angular 
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Figure 13. Drive Control Loop Gain 
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rate information carried as amplitude modulations of the basic oscillator fre- 
quency, and is phase sensitive to the direction of the applied angular rate, the 
drive voltage which was shown to be either in-phase and out-of-phase with 
the angular rate output can be used as the demodulation reference. This is 
shown functionally in Figure 15. The sensor transfer function without damping 
can therefore be represented as 


E q (s) 
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where E q is the transfer function of the demodulated sensor voltage output, 

K q is the sensor scale factor in volts/rad/sec and includes both mechanical 
and electrical gain, co is the separation frequency (difference between vibra- 
tional and output frequencies) and Q is the output quality factor. 


The damping loop, which is also shown functionally in Figure 15, uses the 
damping transducer to apply an amplitude -modulated voltage to the beam- 
transducer to increase the basic damping ratio of the above transfer function. 
The carrier frequency of the damping signal is the basic oscillator frequency. 
The response of the beam- transducers in the output axis due to a damping 
signal at the carrier frequency is 
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This gain must be known to establish the damping loop gain and the effective 
sensor damping ratio. 


The prime design considerations for the damping loop are minimum noise 
and effective control of damping ratio. The passive filter is designed to be a 
minimum bandwidth lead- lag such that sufficient phase lead is obtained at the 
separation frequency. Typically, 
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The modulator is used to modulate the damping error signal in phase with the 
drive voltage and consists of a single transistor chopper . Note that the ac 
error signal produced by the modulator is a square wave. Band- limiting and 
noise suppression is also supplied by the damping amplifier to reduce chopper 
noise. Typically, 


K(s) = 



(107) 


Where co c is the drive frequency and is the midband gain. The sensor 

damping ratio can be adjusted from approximately 0.5 to 1.0 by adjusting the 
damping amplifier gain Kp from 1 to 10. 

Any part of the sensor design which affects the sensor separation or natural 
frequency must be compensated for because of the potential cross-coupling 
and sensor null shifts. The laboratory model was designed using the 
gimballed-flexural pivot approach. This design is such that the inertia of the 
gimbal affects the resonant frequency in one axis but not the other. 

The gimbal inertia acts principally in the readout plane of the sensor and 
tends to lower its resonant frequency with respect to the drive plane, thus 
widening the frequency separation and sensor bandwidth. This is shown in 
the following. 

The first mode resonant frequency for a perfect free-free is given by: 



where: L is the length of the beam 

p is the beam mass per unit length 
El is the flexural rigidity 

and the beam deflection curve is approximately: 

y(x) = Y 0 (sin ^ - f ) ( see Figure 9 ) 
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the beam at the nodes is given by: 


of beam potential and kinetic energies are respectively: 
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Using the above, the resonant frequency can be expressed as: 
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where K„ is the effective rotational spring rate of the beam and I_. is the 
V B 

effective beam rotational inertia. If the gimbal inertia is added at the nodes, 

then the new resonant frequency is: 


*2 

go 



(111) 


where I is the gimbal inertia. 

U 

A sketch of the gimbal is shown in Figure 9 with appropriate nomenclature. 
The inertia is given approximately by: 

Iq - cat ja(a 2 + c 2 ) + 3a 2 + c 2 

where p is density of the gimbal material. 


As an example, the parameters for the design used for preliminary evaluation 
where: 


a = 0.186", c = 0.125", t = 0.08" 
p = 0.3 lb/in 2 , g = 387 in/sec 2 

and 

1^ = 4. 9 x 10 ^ in-lb-sec 2 . 
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For the 4. 0" beam with 0. 172" cross-section, 
Ig = 2. 48 x 10 ^ in-lb-sec^ 

= 4. 5 x 10^ in-lb-rad 
Therefore, 


w = 2140 Hz 


and 


o>* 



2100 Hz 


( 112 ) 


These results are very close to the measured frequency separation which is 
the combination of non-square beam cross-section and gimbal effects. 
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Section 4 


LABORATORY MODEL DESCRIPTION AND TESTING 
DESCRIPTION 


The laboratory model of the free -free beam vibratory angular rate sensor 
was designed with flexural pivot gimballing of the beam and a low-noise set 
of electronics to allow evaluation of the minimum threshold and resolution. 
Pictures of the sensor model are shown in Figures 16 and 17, which show 
an overall view of the beam and electronics, and a close-up of the beam- 
mount to show details. Figure 18 is a sketch of the cross-section of the 
beam and mounting channel at one of the mounting points which shows the 
flexural pivot arrangement. 

The gimballed-beam sensor was selected for fabrication because it repre- 
sented the mounting concept with the best potential for achieving the thres- 
hold and null stability requirements. The combination of the "hard" con- 
struction and the fact that the configuration could be precisely balanced to 
reduce external vibration effects make it superior to either a "soft" or 
rubber -mounted design, or a rigid "clamped -clamped version. The sketches 
of Figures 19 thru 24 give the details and dimensions for the design. The 
pertinent parameters that were selected for the laboratory model were: 


Beam: 

Material 
Length L 
Cross-Section 
Mount Location 

NiSpan C 
= 4. 00 inches 
h y =h z = 0. 172 inches 
x n * = 0. 92 inches 

Gimbal: 

Material 

303 SS 


Cross-Section 

a = 0. 377 inches 


Depth 

C = 0. 172 inches 


Thickness 

t = 0. 082 inches 

Transducers: 

Material 

PZT-4 


Length 

& = 0. 75 inches 


Width 

W = 0. 09 inches 


Thickness 

t = 0. 02 inches 


These parameters were selected to yield beam natural frequencies of ap- 
proximately 2000 Hz and a separation frequency of 50 Hz. Note that the beam 
is made square and that the inertia of gimbal gives the desired separation. 

The mounting locations shown on the sketches are not exact because additional 
length has been left at the beam ends to allow for trimming during test. 
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Figure 16. Solid-State Sensor, Laboratory Model 


Figure 17. Solid-State Sensor Mechanical Assembly 
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TRANSDUCER 


Figure 18. Gimballed Sensor Axes and Transducer Nomenclature * 
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Figure 19. Beam 
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Figure 23. Sensor Box 
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Assembly 


The schematics of the sensor laboratory model electronics are shown in 
Figures 25 thru 28. These electronics are mechanized with 12 microelec- 
tronic operational amplifiers and have the advantage that they are easily 
adaptable to parameter changes during sensor testing. This mechanization, 
although a hybrid with discrete components for setting gains and compensa- 
tion, could be extrapolated to an all microelectronic design. The electronics 
are arranged on four functional boards. 

Figure 25 is the basic drive loop circuit which together with the beam forms 
an amplitude -controlled oscillator that produces a controlled beam vibration 
at its fundamental frequency. The key element is the AGC amplifier which 
is controlled by an FET. Figure 26 consists of the circuits which generate 
the demodulator reference and provide the regulator voltage for the drive 
loop. Figure 27 shows the buffer amplifiers and preamplifier for the sensor 
rate output voltage. The use of dual FET's with a microelectronic amplifier 
reduces the electronic noise present in the output. A quadrature null ad- 
justment is shown before the preamplifier. Figure 28 consists of the demodu- 
lator and damping circuits. A dc null adjustment is included before the output 
amplifier. 

The adjustment and calibration of the electronics for the specific beam being 
used included the following: 

1. Selection of vibration amplitude and adjustment of drive loop gain 
at that amplitude to operate AGC amplifier at approximately the 
midpoint of its gain range. 

2. Quadrature null adjustment after buffer amplifiers. 

3. In-phase null adjustment with drive potentiometer and/or output 
dc null pot. 

4. Adjustment of demodulator reference phase for maximum quadra- 
ture rejection. 

5. Selection of preamplifier gain for desired sensor scale factor. 
LABORATORY MODEL TESTS AND EVALUATION 


The test and evaluation of the laboratory model consisted of the following: 

1. Functional tests of the beam-transducer combination. 

2. Calibration and linearity. 

3. Threshold and resolution measurements. 
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LAYOUT 



Figure 25. Solid-State Sensor Electronics -Board No. 
Drive and Feedback 
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Figure 26. Solid-State Sensor Electronics -Board No. 

Regulation Regulator and Demod. Ref. 




LAYOUT 



Figure 27. Solid-State Sensor Electronics -Board No. 
Readout Buffers and Preamplifier 










4. Drift after start-up and long term drift. 

5. Null drift with temperature. 

6. Frequency response. 

7. Null shift due to vibration. 

8. Cross-coupling tests. 

The functional tests of the mechanical assembly are used for determining 
whether the beam-transducer configuration is of adequate quality and to 
determine accurately various parameters such as resonant frequencies and 
bandwidths. The beam-transducer frequency responses in the drive plane 
and the readout plane are shown in Figures 29 and 30 . The drive or 
vibrational axis for the gimballed beam is the inner axis; the readout axis 
is the outer axis. 

Tests were conducted by exciting the drive (damping) transducer with an 
audio oscillator and monitoring the open or short circuit response of the 
feedback (readout) transducer as a function of excitation frequency. These 
tests are significant because the data provides information about sensor 
quality (transducer gain and bond) proximity of mounting locations to actual 
vibration nodes, and sensor bandwidth and scale factor. Figure 29 is a 
plot of the "Q" curve for the drive plane or axis. The beam gain shown is 
the ratio of the feedback voltage and the drive voltage which are both sinu- 
soidal. The 2KJ2 load is used to simulate short-circuit transducer operation. 
As shown, the "Q” of this unit is unusually high (~ 1300) which is due to the 
accurate location of the beam nodes • 

A similar plot for the readout plane is shown in Figure 30 and indicates 
similar results although the Q is lower (~ 860) due to the extra gimbal and 
pivot effects. The separation frequency or sensor bandwidth for this unit 
is 63. 7 Hz and the scale factor should be between 1. 75 and 2. 00 millivolts 
per degree per second at the readout transducer. 

The data shown in Figures 31 and 32 exhibits the nonlinearities of the beam 
vibrational motion due to the flex-pivot mounting. Since the mounting points 
are fixed rigidly in the frame, large motions of the beam tend to bend the 
pivots laterally which reduces the actual bending motions of the beam and the 
equivalent natural frequency. Figure 31 shows both the amplitude and fre- 
quency nonlinearities. The beam gain changes approximately 2 to 1 over the 
range of drive voltages and the maximum change in resonant frequency is 
approximately 0. 1%. Figure 32 shows the reduction in equivalent Q at 
resonance due to increased beam motion. There is a 50% decrease over the 
range shown. These nonlinearities do not affect the linearity of the unit when 
operated as an angular rate sensor because the drive plane is always 
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- DRIVE VOLTAGE = 10 VOLTS P-P 
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Figure 29. Drive Plane Frequency Response 



DAMPING VOLTAGE = 10 VOLTS P-P 







controlled at a specific operating voltage and the readout plane operates at 
extremely low amplitudes. However, the data shown does imply that the 
long term stability of the sensor will depend upon the changes in the above 
curves with environment and time. Further investigation of these effects 
is planned. 

A summary of the laboratory model parameters is given below: 

Drive plane (inner axis) at 10 volt excitation 


Resonant frequency 2057. 6 Hz 

Q 1290 

Gain (into 2K load) at resonance 0. 08 v/v 

Readout plane (outer axis) at 10 volt excitation 

Resonant frequency 1993. 9 Hz 

Q 865 

Gain (2K load) at resonance 0. 05 v/v 


Separation frequency = Sensor bandwidth = 63. 7 Hz 

The free-free beam sensor was calibrated and its linearity measured for 
input angular rates from 2 deg/sec to 30 deg/sec. The data taken is shown 
in Table 2 and was obtained with a Genesco Rate-of-Turn table with strobe. 
The nominal scale factor was set to be approximately 0.4 volt/deg/sec by 
adjusting the scale factor amplifier gain. For this particular sensor the 
electronic gain was 300 mv dc/mv p-pac and the actual scale factor at the 
sensor was approximately 1. 3 mv p-pac/deg/sec. The data of Table 1 is 
plotted in Figure 33 and indicates an average scale factor of 0. 395/volt/deg/ 
sec. This scale factor was used to compute the point-for -point linearity of 
the sensor and electronics. The linearity data is plotted in Figure 34 and 
represents the point deviation from the average scale factor in deg/sec. The 
data indicates that the linearity is always better than 1%. The sources of 
indicated nonlinearity are null shifts and electronic noise; the beam vibra- 
tions for input angular rates less than 1000 deg/sec are small enough to 
yield better than 0. 1% linearity. 

Some of the most interesting and significant data for the sensor was taken in 
a series of cross-coupling tests which evaluated the null sensitivity of the 
sensor. The objective of this evaluation was to find out to what degree the 
threshold and resolution depended upon the direction of the beam at null. The 
tests were conducted by modulating the drive voltage, and thus the drive 
plane motion, at a frequency below or near the separation frequency and 
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Figure 33. Sensor Output Voltage vs Applied Angular Rate 
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Figure 34. Sensor Linearity 


Angular Rate Sensor Output 

deg/ sec volts 

30 11.16 

25 9.862 

20 7. 906 

15 5.944 

CW 10 3. 958 

8 3. 178 

6 2.380 

4 1. 588 

2 .793 

0 .001 

2 .784 

4 1. 576 

6 2. 374 

8 3. 160 

CCW 10 3. 940 

15 5.910 

20 7.884 

25 9.830 

30 11.05 

Scale factor 0.395 volt/deg/sec 


TABLE 2 
Sensor Calibration 
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monitoring the amplitude of this frequency component in the sensor output 
voltage. In general, the results are a large dependence between sensor 
noise and mechanical null setting. Reductions of up to 10 to 1 in sensor out- 
put noise can be achieved by choosing the correct operating point. Two 
distinct operating points, which are of interest, are the electrical null where 
the sensor output is trimmed using the split drive transducer to 0 volt dc, 
and the mechanical, which is defined as the setting of the two drive voltages, 
which minimizes the noise in the sensor output. Some typical data for the 
sensor is plotted in Figure 35. 

It shows the sensor output voltage due to a modulation of the drive plane 
motion at different settings of the ratio of voltages on the split drive trans- 
ducer. Both axes are calibrated in equivalent deg/sec with an S. F. of 
0. 4 volt/deg/sec. The ratio of sensor coupling at the cross -coupling null 
and at electrical null (0 vdc) is approximately -52 db. Operation of the sen- 
sor at an offset of ~ 16 deg/sec reduces the sensor output noise at the sepa- 
ration frequency by 100 to 1. The separation frequency was chosen for the 
data because the sensor is most sensitive to noise at this frequency. The 
existence of the cross-coupling null is probably due to the high stresses 
placed on the readout crystal due to the drive motion and results in direct 
transmission from the drive transducer to the readout transducer. Further 
work is needed in this area to determine the exact causes of coupling. 

Traces of the sensor output voltage, which are an indication of the sensor 
threshold and resolution, were taken at the various null conditions and for 
just the electronics. These are shown in Figure 36 and show the improve- 
ment obtained by operating at mechanical null. 

Frequency responses for the laboratory model were taken with an oscillating 
rate table for sinusoidal input angular rates of 7 deg/sec peak for the various 
value of sensor damping. This data is shown in Figure 37 and indicates a 
natural frequency of 65 Hz and an inherent damping ratio of 0. 02, which 
agrees very closely with the analytical prediction contained in the Theory of 
Operation. The damping ratio is adjusted by varying the drive amplifier 
gain and, as indicated, a damping ratio of 0. 6 can be achieved. For damping 
ratio of 0. 8 or greater, the larger gains and the chopper mechanization of 
the damping signal increase the level of sensor output noise. 

The laboratory model null characteristics were measured to determine the 
threshold/resolution and the null uncertainty over a long period. This type 
of data is shown in Figures 38 and 39. The short-term random fluctuations 
of the sensor output voltage, similar to those exhibited during the discussion 
of the cross-coupling characteristics, define the sensor threshold or signal- 
to-noise ratio. Since the vibrating beam sensor does not have any type of 
deadband caused by friction or other causes, the threshold depends on 
noticeable changes in the output noise. The trace in Figure 38 indicates a 
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FEEDBACK 10 VOLTS P-P 
DRIVE 4 VOLTS P-P 

DRIVE MODULATION 0.4 VOLTS P-P 
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Figure 36. Sensor Null Traces 
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Figure 39. Null Drift of Sensor 


peak-to-peak noise voltage of approximately 0. 7 millivolts which with a sen- 
sor scale factor of 0. 4 volt/deg/sec is equivalent to ~ 0. 002-0. 003 deg/sec. 
The primary source of the random fluctuations is the electronic noise which 
results because of the large amplifier gains required. Another source is 
random modulations of the drive motion which produces separation frequency 
noise at the output; this source has been minimized by operating at mechanical 
null. The long term type of data shown in Figure 39 is the null shift of the 
laboratory model after turn-on. As shown, the sensor has an hourly type of 
exponential drift with an approximate time constant of 4 to 5 hours. The 
amplitude of this drift is about 0. 15 deg/sec. The cause of this drift is not 
known exactly but is probably due to a stabilization of transducer and mount 
parameters. This drift is a function of the time that the sensor has been 
turned on. 

The environmental testing of the sensor laboratory model consisted of tem- 
perature tests of the mechanical assembly from -10°C to +70°C and vibration 
tests of the sensor at lg from 10 to 2000 Hz. The results of the temperature 
tests are shown in Figure 40 which is a plot of sensor null in equivalent 
deg/sec versus temperature. For these tests, the beam and mounting struc- 
ture were inside the temperature oven and the electronics were outside. The 
readings at each temperature were made after a half hour soak. As indicated 
in Figure 40, the null shift due to temperature is approximately ± 6 deg/ sec 
for -10°C to +70°C. The principal cause of these shifts is thought to be the 
drive transducer and its interface. As pointed out in the section on Theory 
of Operation, the sensor null is particularly sensitive to changes in the 
homogeniety of drive transducer coupling. The data shown in Figure 40 con- 
sists of a typical run from -10°C to +70 or C and the band of data points from 
five such runs. The variation from run to run is probably due to changes in 
the transducer characteristics. In addition to this data, the scale factor 
was checked at various temperatures and was found to vary approximately 
± 10% over this range. 

The vibration data is shown in Figure 41. The mechanical assembly was 
vibrated in the output axis which is the most sensitive axis of the sensor to 
vibration inputs. At low frequencies, below beam resonance, the null shift 
in equivalent deg/sec due to lg vibration is less than 0. 25 deg/sec. At beam 
resonance, the null shift increases at least an order of magnitude or two due 
to Q multiplication and indicates that the nodal points were not exactly cor- 
rect for this beam. The accuracy of the data at or near resonance is question- 
able because of the frequency tolerances of the vibration input. 
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Figure 41. Sensor Vibration Response 


Section 5 


CONCLUSIONS AND RECOMMENDATIONS 


The development and design of the laboratory model for the "free-free" beam, 
as well as the analysis and testing done on this contract, has shown speci- 
fically that improvements in the design of the beam mounting structure and 
the electronics can improve the null characteristics of the sensor by reducing 
noise output and stabilizing it. This is accomplished, primarily, by elec- 
tronically decoupling the vibrations in the sensor and mechanically decoupling 
the beam vibrations from external accelerations. These improvements are 
the first step toward the development of a vibratory angular rate sensor of 
inertial quality. The work on a complete mathematical model for the free- 
free beam sensor is a start toward the optimization of all sensor design 
parameters for minimum threshold and null uncertainty. 

The "free-free" beam vibratory angular rate angular rate sensor in its pre- 
sent development status is applicable to a wide range of control systems 
requiring angular rate measurement. These include the damping sensor for 
a mass expulsion attitude control system for spacecraft and the stability 
augmentation sensor for aircraft flight control. In particular, for applications 
where constant rate measurement is not needed and a "canceller" or high 
pass filter can be used, the sensor can be considered for immediate applica- 
tion. The filter, which may have a time constant of 60 seconds or larger, 
serves to eliminate the null uncertainty and drift due to turn-on or tempera- 
ture changes. 

Further development of the sensor to attain lower thresholds and null un- 
certainties should include an investigation of the use of magnetic drive and 
readout of the beam vibrations, as well as material research. The develop- 
ment of the flexural pivot mount has made it feasible to achieve dimensional 
stability between the beam and its mounting structure, and thus has made a 
magnetic drive possible. The principle improvement in sensor performance 
to be gained from a magnetic mechanization is long-term null stability. 
Continued development should also include a parameter optimization study 
and further work on a complete math model for the sensor. The work done 
on this contract has shown the interrelationships between the sensor design 
parameters, and further work is needed to be able to systematically select 
dimensions, frequencies, and gains for minimum sensor threshold. 


84 



Appendix A 

FREE LATERAL VIBRATION OF A FREE-FREE BEAM 

The partial differential equation for the deflection x(z, t) of a vibrating beam 
of uniform cross-section such as shown in Figure A- 1 is 



Figure A-l. Free-Free Beam 

EI ? — _? — ^(M) (page 324 of reference 1) 

*z 4 dS 

where j u is mass of the beam per unit length and EI is flexural rigidity. The 
solution of the above equation is found by assuming a solution of the form 

x(z, t) = H(z) sin cot 

Therefore 
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or 


4 

EI — sin = +co 2 p. H(z) sin cot 

dz 4 


Define 

then 


The general solution of the preceding equation is found by Laplace transform 
techniques and is 

H(z) = A t cos kz + Ag sin kz + cos h kz + A^ sin h kz 
where the A's are arbitrary constants. 

The solution of this problem for the particular case of the free -free beam is 
found by applying the boundary conditions 

= 0 
z = 0 

= 0 
z = 0 

= 0 
z = St 

= 0 
z = ft 


(1) 

d 2 H(z) 
dz 2 _ 

(2) 

d 3 H(z) 

. 3 

dz _ 

(3) 

d 2 H(z) 
_ dz 2 _ 

(4) 

d 3 H(z) 

. 3 
dz 


d 2 H(z) 

dz 4 

, 2 

k 4 - 

- EI 


2 

co ii 

EI 


H(z) = 0 


5^1 - k 4 H(z) = 0 

dz 
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The conditions above state that the shape of the beam at any time t must be 
such that the bending moment and shearing force vanish at the ends. By 
applying conditions (1) and (2), we find that 

Aj = Ag and A ^ 


so that 


H(z) = Aj (cos kz + cos h kz) + Ag (sin kz + sin h kz) 


Applying conditions (3) and (4), we obtain the two equations 

Aj (cosh kf - cos kf ) + Ag (sinh kf - sin kf ) = 0 

Aj (sinh kf + sin kf ) + Ag (cosh kf - cos kf ) = 0 


A non-trivial solution of these equations for the constants A^ and Ag requires 
that the determinant of the coefficient matrix be equal to zero. 

(cosh kf - cos kf ) 2 - (sinh 2 kf - sin 2 kf ) = 0 
or cos kf cosh kf = 1 


The roots of this equation are the eigenvalues of the problem which deter- 
mine the resonant frequencies of the beam. The first six roots of the 
equation are 


For 


For 


kf = 0, 4.730, 7.853, 10.996, 14.137, 17.279 
kf = 0, k = 0 and u> = k 2 j 1/2 = 0 

H(z) = 2Aj = constant 


x(z, t) = H(z) sin cot = 0 


4 73 

kf = 4.73, k = , co 


22.4 




1/2 


^2 _ cosh 4. 73 - cos 4.73 
Aj ~ sin 4.73 - sinh 4.73 


87 


TW x . , 4. 73z _ L 4. 73z 

H(z) = (cos — j — + cosh — j — 


- sin 


4. 73z 

£ 


- sinh 


4. 73z 

£ 


) 


The above shape of the beam for the first bending mode as given by H(z) is 
the normal elastic curve for a free-free beam. Note that the curve has a 
translation and rotation with respect to the z-axis due to constant and first 
degree terms in z which are also solutions of the partial differential equation. 

An approximation to the elastic curve with the rotational terms removed is 

H(z) = x sin -r- - a 
o £ 

The approximate solution for x(z,t) for the first mode is 
x(z,t) = H(z) sin cut 
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x(z,t) = (x Q sin — - a) sin cut 


where 


co = 


22.4 



1/2 


and a is a constant that determines the location of the nodes. 


The actual value of ’a' can be found from the fact that since no external 
alternating force is acting on the beam, its total vertical momentum must 
be zero. When the beam passes through its equilibrium position, the ends 
have downward velocities and the middle has an upward velocity. The total 
momentum is zero for all t if the areas above and below the z-axis are equal. 


For 


cot = 


£ 

2 


0 = 



x(z)dz = x 

o 




Therefore, 


a 



and 


x(z,t) =x Q 




sin cot 
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The nodes are located at the points where 


or z = 0. 22* 

z = 0.78* 
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